main scVelo steps

first issue, cell barcode names to match

all kept

stochastic model

seems somehow reasonable, Ccl9+ circulating-like EOS3 to=> Nmur1- EOS2 to=> Nmur1+ EOS3,
and the spliced-unspliced ratio supports,
but local stream looks like a little strange that does not totally fits the final prediction

dynamic model

same final prediction but similar chaotic local stream.
this analysis is only for a quick test,
further modification may need more detailed demonstrations on gene level and ranges of cutoffs
to test the robustness and differences